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A 2.25 x 10 s J/tp event sample accumulated with the BESIII detector is used to study the doubly 
OZI suppressed decay modes J/tp —t "fu<j>, w — > ir + ir~ir°, <p — > K + K~ . A strong deviation (> 30<r) 
from three-body J/tp — ► -yuicp phase space is observed near the u}(j> mass threshold that is consistent 
with a previous observation reported by the BESII experiment. A partial wave analysis with a tensor 
covariant amplitude that assumes that the enhancement is due to the presence of a resonance, the 
X(1810), is performed, and confirms that the spin- parity of the X(1810) is ++ . The mass and width 
of the X (1810) are determined to be M = 1795 ± 7(stat)t^ 3 (syst)±19(mod) MeV/c 2 and V = 95 ± 
10(stat) ^34(syst)±75(mod) MeV/c 2 , respectively, and the product branching fraction is measured 
to be B(J/tp -> 7^(1810)) x B(X(1810) u<j>) = (2.00 ± 0.08(stat)±?;^(syst)±1.30(mod)) x 1(T 4 . 
These results are consistent within errors with those of the BESII experiment. 

PACS numbers: 14.40. Lb, 14.40.Be, 13.25.Gv 

I. INTRODUCTION 

An anomalous near-threshold enhancement in the ivcp invariant-mass spectrum in the process J/ijj — > juxp w &s 
reported by the BESII experiment [1]. A partial wave analysis (PWA) that used a helicity covariant amplitude 
that assumed that the enhancement was produced by a resonance, denoted as the A"(1810), was performed on the 
BESII event sample. The analysis indicated that the A"(1810) quantum number assignment favored J PC = ++ 
over J PC = h or 2 ++ with a significance of more than 10a. The mass and width were determined to be M = 
1812±2§(stat)±18(syst) MeV/c 2 and T = 105 ± 20(stat)±28(syst) MeV/c 2 , respectively, and the product branching 
fraction B(J/tp 7 X(1810)) • B(X(1810)-+ uj<I>) =[2.61 ± 0.27(stat)±0.65(syst)] x 10" 4 was measured. The J/tp -> 
"fW(p decay mode is a doubly OZI suppressed process with a production rate that is expected to be suppressed relative 
to J/tp — > juLd or J/tp —> 7</></> by at least one order of magnitude [2]. Possible interpretations of the ojcp threshold 
enhancement include a new type of resonance, such as a tctraquark state (with structure q 2 q 2 ) [3], a hybrid [4], or 
a glucball state [5] etc., a dynamical effect arising from intermediate meson rescattcring [6], or a threshold cusp of 
an attracting resonance [7]. As of now none of these interpretations has either been established or ruled out by 
experiment. 

A search for the X(1810) was performed by the Belle collaboration in the decay of B^ 1 — > K^ujcp [8], but no 
obvious AT(1810) signal was observed. A high statistics data sample collected with the BESIII detector provides a 
good opportunity to confirm the existence of the u<p threshold enhancement, study its properties and search for other 
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possible related states that decay to u<j>. 

In this paper we present a PWA that uses a tensor covariant amplitude for the J/ip — > 7W0 process, where the 4> 
is reconstructed from K + K~ and the w from 7r + 7r~7r°. The analysis is based on a sample of (225.3 ± 2.8) x 10 6 J/ip 
events [9] accumulated with the new Beijing Spectrometer (BESIII) [10] located at the Beijing Electron-Positron 
Collider (BEPCH) [11]. 

II. DETECTOR SETUP AND MONTE CARLO SIMULATION 

BEPCII is a double-ring e + e~ collider designed to provide a peak luminosity of 10 33 cm' 2 s _1 with beam currents 
of 0.93 A. The BESIII detector has a geometrical acceptance of 93% of 4-7T and has four main components: (1) A 
small-cell, helium-based (40% He, 60% C3H8) main drift chamber (MDC) with 43 layers providing an average single- 
hit resolution of 135 /im, charged-particle momentum resolution in a 1 T magnetic field of 0.5% at 1 GeV/c, and 
a dE/dx resolution better than 6%. (2) An electromagnetic calorimeter (EMC) consisting of 6240 CsI(Tl) crystals 
in a cylindrical structure (barrel) and two endcaps. The energy resolution for 1.0 GcV/c 7-rays is 2.5% (5%) in the 
barrel (endcaps), and the position resolution is 6 mm (9 mm) in the barrel (endcaps). (3) A time-of-flight system 
(TOF) constructed of 5 cm thick plastic scintillators, with 176 detectors of 2.4 m length in two layers in the barrel 
and 96 fan-shaped detectors in the endcaps. The barrel (endcap) time resolution of 80 ps (110 ps) provides 2a K/ir 
separation for momenta up to ~ 1.0 GeV/c. (4) The muon system (MUC) consists of 1000 m 2 of Resistive Plate 
Chambers (RPCs) in nine barrel and eight endcap layers and provides 2 cm position resolution. 

In this analysis, a GEANT4-based [12] Monte Carlo (MC) simulation software package, BOOST [13], is used. It 
provides an event generator, contains the detector geometry description, and simulates the detector response and signal 
digitization. The production of the J/i/j resonance is simulated by the Monte Carlo event generator KKMC [14, 15], 
while the decays are generated by BesEvtGen [16, 17] for known decay modes with branching ratios set at the PDG [18] 
world average values, and by the Lund-Charm model [19] for the remaining unknown decays. The analysis is performed 
in the framework of the BESIII Offline Software System (BOSS), which takes care of the detector calibration, event 
reconstruction, and data storage. 
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III. EVENT SELECTION 

Signal J/ip —> juxj) events with 4> — > K + K~ and lu —> iT + ir~ir final states have the topology 3jK + K~tt + it~ . The 
event candidates are required to have four well reconstructed charged tracks with net charge zero, and at least three 
photons. 

Charged-particle tracks in the polar angle range |cos6*| < 0.93 are reconstructed from the MDC hits. Only the 
tracks with points of closest approach to the beamline that are within ±10 cm of the interaction point in the beam 
direction, and within 1 cm in the plane perpendicular to the beam are selected. TOF and dE/dx information are 
combined to form particle identification confidence levels for 7r, K and p hypotheses. Kaons are identified by requiring 
the particle identification probability (Prob) to be Prob(K) > Prob(Tr) and Prob(K) > Prob(p). Two identified 
kaons with opposite charges are required. 

Photon candidates are reconstructed by clustering signals in EMC crystals. The energy deposited in the nearby 
TOF counters is included to improve the photon reconstruction efficiency and its energy resolution. The photon 
candidates are required to be in the barrel region (|cos#| < 0.8) of the EMC with at least 25 MeV total energy 
deposition, or in the endcap regions (0.86 < | cos#| < 0.92) with at least 50 MeV total energy deposition, where 9 
is the polar angle of the shower. The photon candidates are furthermore required to be isolated from all charged 
tracks by an angle > 10° to suppress showers generated by charged particles. The showers in the region between the 
barrel and the endcaps of the EMC are poorly measured and excluded. Timing information from the EMC is used to 
suppress electronic noise and energy deposits that are unrelated to the event. Events with at least three good photon 
candidates are selected. 

A four-constraint (4C) energy-momentum conserving kinematic fit is performed to the 3~fK + K~Tr + 7r~ hypothesis. 
For events with more than three photon candidates, the candidate combination with the minimum x\c 1S selected, 
and it is required that x\c < 40 (the requirement is determined by optimizing S/y/ S + B, where S is the number of 
MC signal events generated with phase space, and (S + B) is the number of signal plus background candidate events 
in the data). In order to remove background stemming from J/ip — >2jK + K~tt + tt~ and J/ip — >A"fK + K~tt + tt~ , we 
performed 4C kinematic fits for the hypotheses of 2jK + K ~ir + ir~ and ^K + K ~ir + ir~ (for the events that have at least 
four good photon candidates). We require x1c(3"fK + K~'K + ir~) < xlc (^jK + K^tt + tv^) and XicfilK + K~7r + ir~) < 
x\q{^K + K~ir + ir~), respectively. The 7r° candidates are reconstructed from the two of the three selected photons 
with invariant mass closest to the ir° mass, and |M 77 — M n o\ < 20 McV/c 2 is required. 
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FIG. 1: (a) A scatter plot of M K + K - versus M^+^-^o. The boxes indicate the signal region labeled as S and sideband regions 
labeled as A, B and C (defined in text), (b) The K + K~ invariant-mass distribution; the shaded histogram shows the events 
within the uj sideband region, (c) The ■k + tt~tv° invariant mass distribution; the shaded histogram shows the events within the 
cj> sideband region. 
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FIG. 2: (a) The jty + tt ty° invariant-mass distribution, (b) The K + K tt + ty 7t° invariant-mass distribution; the dashed 
line is the mass distribution of the phase space MC sample; the solid histogram shows the mass distribution without the 
M(77T + 7r _ 7r°)>I.0GeV/c 2 requirement, (c) A Dalitz plot of M 2 (77r + 7r~7r ) versus M 2 (~{K + K~). 



A scatter plot of the M^+k- versus M w +„-„o invariant masses for events that survive the above selection criteria 
is shown in Fig. 1(a). One cluster of events populates the <p4> region, which arises from the well known process 
J/ip — > ~f4>4> (one 4> — > 7T + 7r~7r°, the other <f> —> K + K~), and another cluster of events shows up in the ui(f> signal 
region. Since the decays of J /ip — > cjcA and J/tp — > tt ujcj) are forbidden by C parity conservation, the observed 
events in the uj(f> region are an unambiguous signal for the radiative decay process J/ip — > ^u<j). The mass window 



requirements (T)\M„+^-„o — M u \ < 40McV/c 2 (the requirement is determined by optimizing S/y/S + B) and (II) 
\M K + K - — MA < 15McV/c 2 (the requirement is determined by optimizing S/\/S + B) arc defined for the w and 
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signal region, respectively, while the requirements of (III) 60McV/c 2 < \M n + n - n o — M u \ < 140MeV/c 2 and (IV) 
1045MeV/c 2 < M K + K - < 1075 MeV/c 2 are defined for the ui and sideband regions, respectively. Figure 1(b) shows 
the K + K~ invariant-mass distribution for events in which the 7r + 7r _ 7r° invariant-mass lies within the ui signal range 
(requirement I); here a signal can clearly be seen. The shaded histogram in Fig. 1(b) shows the corresponding 
distribution for events within the u> sideband region (requirement III). A small signal from the J /ip —> 707r + 7r~7r° 
background is evident. Figure 1(c) shows the tt + tt~tt invariant-mass distribution for events with K + K~ invariant- 
mass within the signal range (requirement II). As expected, uj and signals are clearly seen. A small r\ signal is 
also observed; this comes from the decay chain J/ip — ¥ r yqK + K~ (rj — > 7r + 7r _ 7r°). The shaded histogram in Fig. 1(c) 
shows the corresponding distribution for the events within the sideband region (requirement IV). For events that 
survive the to and requirements on the 7r + 7r _ 7r° and K + K~ invariant mass (requirements I and II), respectively, 
the 77r + 7T~7r invariant- mass distribution is shown in Fig. 2(a). Here an r) 1 peak is observed; this comes from the 
decay chain J/ip — > 07/(0 — > K + K~ , rf — > jco, uj — > 7t + it~it ). To characterize these events, a large MC sample of 
J/tp — > 4>ri' is generated with a flat angular distribution. These have a K + K~ir + Tt~ir® invariant mass distribution 
that is concentrated at masses higher than 2.5GeV/c 2 and have no impact on the uj(j> mass threshold region of interest. 
A further requirement M(77r + 7r _ 7r°)>1.0GcV/c 2 (requirement V) is imposed to remove background from J/ifj — > <j>rf . 
Figure 2(b) shows the invariant mass of K + K~TT + Tt~TT {) for events with requirements I, II and V applied, where a 
peaking structure near the u>4> invariant-mass threshold is observed. The solid histogram in the figure shows the 
K + K~ir + ir~ir° invariant-mass distribution without requirement V. The invariant mass distribution is very different 
from a pure phase-space distribution from MC (dashed histogram, arbitrarily scaled) . The threshold structure shows 
up as a diagonal band along the upper right-hand edge of the Dalitz plot in Fig. 2(c). 

The observed us<f> mass-threshold enhancement is similar to that observed by the BESII experiment [1]. To ensure 
that the enhancement is not due to some background process, detailed studies of potential background sources have 
been performed using both data and MC. Non-w and non-0 backgrounds are studied using ui and <f> sideband data. 
Figure 3(a) and 3(b) show the K + K ~ir + ir~ir° invariant mass for events in the ui sideband region (labeled as Box 
A in Fig. 1(a)) and the </> sideband region (labeled as Box B in Fig. 1(a)); these are used to determine the non-w 
and non-0 background contamination in the signal regions. Figure 3(c) shows the same distribution for events in the 
corner region (labeled as Box C in Fig. 1(a)), for which both the K + K~ and the 7r + 7r~7r° invariant masses are in 
the 4> and uj sidebands; these are used to estimate the non-0 non-w background. The background contamination in 
the signal region is estimated to be the sum of the Fig. 3(a) and Fig. 3(b) sideband distributions with the Fig. 3(c) 



9 



> 
C5 

S 20\ 

5a 
c 

w 101 



,.1 



2.5 



(a) 



i; 



M(K + KVjfit°)(GeV/c') 



> 
o 



5a 

w 10 



2.5 



(b) 



II J 



M(K + KVir-ir°)(GeV/c 2 ) 




2 2.5 3 2 2.5 3 

M(K + KVifit°)(GeV/c 2 ) M(K + K"jr + jur )(GeV/c 2 ) 



FIG. 3: The K + K - 7r + 7r~7r° invariant-mass distribution for (a) the events in the uj sideband region (Box A in Fig. 1(a)); (b) 
the events in the <j> sideband region (Box B in Fig. 1(a)); (c) the events in the corner region (Box C in Fig. 1(a)); (d) the events 
in the ucj) signal region; the solid histogram is the background distribution estimated from the sideband events, the dashed 
histogram is that obtained from inclusive J/ip MC samples. 

distribution subtracted to account for double counting of non-0 non-w background in Fig. 3(a) and Fig. 3(b). Phasc- 
space-MC-determined normalization factors are applied that account for differences in the sizes of the selected regions 
and the difference in the available phase space in the signal and sideband regions. The background contamination 
in the signal region determined in this way is shown as a solid histogram in Fig. 3(d). The shape of the estimated 
background is very different from that of data in the signal region, and no evidence of an enhancement near the cucf) 
mass threshold is observed from the non-w and non-0 background events in the data. 

An inclusive MC sample of 225M J/ip events generated according to the Lund-Charm model [19] and the PDG decay 
tables is also used to study the potential backgrounds. The dashed histogram in Fig. 3(d) shows the K + K~Tr + ir~Tr a 
invariant-mass distribution for the selected inclusive J/ip MC events, where no peaking background at the ui<j) invariant- 
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mass threshold is observed. Exclusive background MC samples of J/tp decays that have similar final states are 
generated to further investigate possible background sources. The main backgrounds come from J/ip — > loK*K 1 K* — > 
Kn° and J/ip -> w/i(1420), /i(1420) -> K + K~n° events. For these, the K + K ~ir + TT~ir° invariant mass distribution 
peaks at high masses, and none of them channels produce peaking structures at the unp mass threshold. 



IV. PARTIAL WAVE ANALYSIS 



A P WA was performed on the selected J/ip — > ^uj<p candidate events to study the properties of the uxp mass threshold 
enhancement. In the PWA, we assume the enhancement is due to the presence of a resonance, denoted as X, and 
the decay processes are described with sequential 2-body or 3-body decays: J/ip — > 7X, X — > uj<j), 10 — > 7r + 7r~7r° and 
<j> — > K + K~ . The amplitudes of the 2-body or 3-body decays are constructed with a covariant tensor amplitude 
method [20]. The intermediate structure X is parameterized with the Breit-Wigner propagator 

BW = \/{M 2 - s-iMY) (1) 

with constant width, where s is the u<j) invariant mass-squared, and M & T are the resonance mass and width, 
respectively. The amplitude for the sequential decay process is the product of all decay amplitudes together with the 
Breit-Wigner propagator. The total differential cross section da/d& for the process is the square of the linear sum of 
all possible partial wave amplitudes: 

% = \J2A(J™)\\ (2) 

where A(J PC ) is the total amplitude for all possible resonances with given J PC . 

The relative magnitudes and phases of the states arc determined by an unbinned maximum likelihood fit of the 
measured cross section da/d§. The basis of likelihood fitting is the calculation of the probability that a hypothesized 
probability distribution function can produce the data set under consideration. The probability to observe the event 
characterized by the measurement £j is the differential cross section normalized to unity: 

- rftfly (3) 

where £»;(£») = (g§)i and e(£i) is the detection efficiency. The joint probability density for observing the N events in 
the data sample is: 

^n^njfgfo. (4) 
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FUMILI [21] is used to optimize the fit parameters in order to achieve the maximum likelihood value. Technically, 
rather than maximizing C, S = -ln£ is minimized, i.e., 



In practice, the normalized integral J d£iw(£i) e (&) ^ s evaluated using the J/ip — > ^ucj) phase space MC sample. For a 
given data set, the second term is a constant and has no impact on the relative changes of the S value. The details 
of the PWA fit process are described in Rcf. [22]. In the minimization procedure, a change in log likelihood of 0.5 
represents a one standard deviation effect for the one-parameter case and is used to evaluate statistical errors. 

Conservation of J PC , in the J/ip —> jX,X ujcp process in the case of a pseudoscalar intermediate resonance X, 
allows only V wave contributions in both the radiative decay J r /ip — > 7 A and the hadronic decay X — > u<p. For the 
production of a ++ , 1 ++ or 2 ++ resonance, both S and T> waves are possible for both the radiative and hadronic 
decays, but only the S wave contribution is considered in the fit, since the D wave can be expected to be highly 
suppressed near the mass threshold. Intermediate X structures with J = 2 h or higher spin are not considered in 
the analysis. To investigate the J PC of the X(1810), we tried different J PC assignments in the fit, and the assignment 
with the best log likelihood value is identified as the J PC of the X(1810). Some known mesons, e.g. /2(1950) or 
/o(2020), with a mass above the uxp invariant-mass threshold, are expected to decay to w<p final states. All possible 
mesons listed in the PDG tables are included in the fit. To consider the contribution from phase space, J/ip — > ^uxp 
without an intermediate state X is also included in the fit with an amplitude modeled by the same sequential process 
and a very broad intermediate state, i.e., M = 2500MeV and T = 5000MeV. In the PWA fit, the phase space is assigned 
to a given J PC , which is determined by the optimization of the likelihood fit. The background event contribution to 
the log likelihood value is estimated from the weighted events in the sideband region, and subtracted in the fit. 

In the PWA fit, different J PC combinations of the X(1810) structure and the phase-space contribution, as well as 
different combinations of additional mesons listed in the PDG tables, are tried. The mass and width of the A(1810) 
are determined by a scan of the maximum log likelihood value, while the mass and width of the additional mesons are 
fixed with their PDG values. The statistical significance of the state is determined by the changes of the maximum log 
likelihood value and of the number of degrees of freedom (And/) in the PWA fits with or without the state included. 
Only states with statistical significance larger than 5a are included in the best solution. 

Finally, together with the contributions of the A(1810) and phase-space, additional ++ , 2 ++ , and h components 
are found (> 5a) in the best solution of the PWA fit. In the following, the masses and widths of the ++ , 2 ++ and 
h components are assigned to be those of /o(2020), /2(1950) and 7?(2225), respectively, since the fit with these has 




(5) 
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the best log likelihood value. Various PWA fits with different ++ , 2 ++ and h components were also performed. 
The results for the A(1810) are robust, while the fit is not very sensitive to the masses and widths of the ++ , 2 ++ 
and h components. The log likelihood values changes are rather small when the /o(2020), /2(1950) and 77(2225) are 
replaced by other resonances with the same J PC and similar masses. The details are shown below. The J PC = ++ 
assignment for the A(1810) has by far the highest log likelihood value among the different J hypotheses. The 
minus log likelihood value (S) for a J PC = ++ assignment to the (1810) is 227 below that of the second lowest 
value (obtained for a J PC = 2 ++ assignment), and is 783 below that for a fit with the AT(1810) omitted. The latter 
corresponds to a statistical significance of more than 30(7. Different J PC assignments for the phase space contribution 
are tested in the PWA fit and J PC — h is favored. The assigned values for the J PC , mass, width and number of 
events for the five components for the best fit solution arc summarized in Tabic I. The mass and width of the X(1810) 
are obtained to be M = (1795 ± 7) MeV/cr and T = (95 ± 10) MeV/c 2 , respectively, where the errors are statistical 
only. The contributions of each component of the best solution of the PWA fit are shown in Fig. 4(a). The changes of 
the log likelihood value AS and of the number of degrees of freedom And/ that occur when a state is dropped from 
the PWA fit, as well as the corresponding statistical significance, are also listed in Table I. The statistical significance 
of the /2(1950), /o(2020) and ?y(2225) contributions are 20.4cr, 13.9c and 6.4er, respectively. The reconstruction 
and final-selection efficiency of the AT(1810) is determined from a weighted phase space MC sample of J/ip — > 
where the weight is the differential cross section for the measured events calculated with the magnitudes and phases 
of the partial amplitudes from the best solution of the PWA fit. The efficiency is determined to be 6.8% and the 
corresponding branching fraction is B{J/ip -> 7^(1810)) x B(X(1810) -> cj0) = (2.00 ±0.08) x 10~ 4 , where the error 
is statistical only. 



TABLE I: Results from the best PWA fit solution. 



Resonance 


J PC 


M(MeV/c 2 ) 


T(MeV/c 2 ) 


Events 


AS 


And/ 


Significance 


X(1810) 


0++ 


1795 ± 7 


95 ±10 


1319 ± 52 


783 


4 


> 30a 


fa (1950) 


2++ 


1944 


472 


665 ± 40 


211 


2 


20.4cr 


f (2020) 


0++ 


1992 


442 


715 ±45 


100 


2 


13.9o- 


n(2225) 


0"+ 


2226 


185 


70 ±30 


23 


2 


6.4a 


phase space 


o-+ 






319 ± 24 


45 


2 


9.1a 



The invariant-mass spectra M(K + K ir + ir ir°), M(^tt + tt it ), M('jK + K ) and the cos# 7 , cos^, cos9<p, cos9k, 
00, and x angular distributions of the data and the PWA fit projections with the best solution as well as the different 
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components are shown in Fig. 4. Here the angles # 7 , 6 U , 9^ and Ok are the polar angles of the radiative photon in 
the J/ip rest frame, the normal to the to decay plane in the uj system, <f> meson momentum direction in the uxj) rest 
system, and the kaon from <f> decay in the 4> rest system, respectively; 00 is the azimuthal angle of the 4> meson in the 
oj(f> system and x is the angle between azimuthal angles of the normal to the ui decay plane and the momentum of a 
kaon from <f> decay in the uj(f) system. The PWA fit projection is the sum of the signal events with the best solution 
and the background estimated from the weighted events in the sideband region. 

To determine the goodness of fit, a \ 2 is calculated by comparing the data and fit projection histograms, where x 2 
is defined as [22] 

f-f Vi 

and n,; and Vj are the number of events for the data and the fit projections with best solution in the i th bin of each 
figure, respectively. The \ 2 an d the number of degrees of freedom (ndf) for each mass and angular distribution are 
shown in Table II, where the number of bins is taken as the number of degrees of freedom. The values of x 2 l n( M 
range between 0.62 and 1.70, indicating reasonable agreement between data and the fit. 



TABLE II: Goodness of fit check for the invariant-mass distributions and angular distributions shown in Fig. 4 



Variable 




M(77r + 7T~7r°) 


M(jK+K~) 


0- f 


Oui 


0,i, 


9k 


4><t> 


X 


x 2 


AAA 


36.4 


42.4 


24.2 


12.4 


28.2 


18.2 


26.4 


51.0 


ndf 


40 


35 


40 


20 


20 


20 


20 


20 


30 


X 2 /ndf 


1.11 


1.04 


1.06 


1.21 


0.62 


1.41 


0.91 


1.32 


1.70 



Additional fits with different assumptions have been carried out to check the influence on the parameters of the 
X(1810). The masses and widths for the /2(1950), /o(2020) and ^(2225) are difficult to determine accurately from 
this analysis and the achieved accuracy can not compete with the PDG accuracy because of the dominant X (1810) 
component; instead they are fixed to their PDG values in the fits. If we change the masses and widths of these three 
mesons by one standard deviation in the fitting, the log likelihood value changes by AS < 3 after refitting mass and 
width of the A(1810); these values and the branching fraction remain consistent within the statistical errors. The 
maximum difference is taken as a systematic error. An alternative method to test the influence of the parameters of 
the three known mesons is to replace the mesons by states listed in the PDG tables with the same J PC and similar 
mass. When the parameters of /2(1950) are replaced with those of the /2(1910), the log likelihood value increases by 
AS = 6 after refitting the mass and width of the A(1810). When the parameters of /o(2020) are replaced with those 
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FIG. 4: Comparisons between data and PWA fit projections: (a) The K + K~ ir + ir~ tv° invariant-mass distribution; (b) the 
77r + 7T - 7r invariant-mass distribution; (c) the r yK + K~ invariant-mass distribution; (d) the polar angle of the radiative photon 
(0 7 ); (e) the polar angle of the normal to the ui decay plane in the to system (0J)\ (f) the polar angle of the c/> in the uxf) rest 
system (6^,); (g) the polar angle of the kaon in the <f> rest system (6k)', (h) the azimuthal angle of the <j> in the uxf) system; (i) 
the distribution of \ which is the angle between azimuthal angles of the normal to the uj decay plane and the momentum of a 
kaon from c/> decay in the ujcj) system. 
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of the /o(2100), the log likelihood value increases by AS = 13. If the /o(2020) is replaced with ++ phase space, the 
log likelihood value increases by AS = 9. If the /2(1950) is replaced with /2(1910) and the /o(2020) is replaced with 
a ++ phase space, the log likelihood value increases by AS = 10 after refitting. By comparing the log likelihood 
values, the combination of A(1810), /2(1950), /o(2020), 7^(2225) and J PC = h for the phase-space contribution is 
found to be the best solution, and the mass, width and branching fraction of X(1810) changes are less than twice the 
statistical errors. 

The states listed in the PDG tables with mass above the uicj) threshold that are consistent with decaying into uicf) 
under spin-parity constraints, arc the / 2 (1910), / 2 (2010), /o(2100), / 2 (2150), /o(2200), / 2 (2300) and / 2 (2340), etc. 
Relative to the best solution of the PWA fit, as these resonances are added in the fit, the log likelihood value improves 
by 11.7, 8.1, 1.5, 3.2, 1.7, 2.5 and 0.9 after refitting the mass and width of the A(1810), and the statistical significance 
of these additional resonances are all less than 5er, while the mass, width, and branching fraction of X(1810) are 
consistent with those from the best solution within statistical errors. The maximum difference between the best fit 
result and the result with extra states included is taken as a systematic error. In the best fit solution, phase space is 
included and approximated as a broad J = h resonance. An additional phase-space distribution amplitude with 
different J PC was added to test whether the data contain different J PC phase-space contributions. When the fit is 
redone including additional phase-space contributions with J PC = ++ , 1 ++ , 2 ++ , the log likelihood value improves 
by 0.1, 3.8 and 3.7, respectively. No evidence of phase space contributions with different J PC values is found, while 
the X(1810) mass, width, and branching fraction are consistent with best solution within statistical errors. The 
maximum differences are taken as systematic errors. The BESIII collaboration has observed two new pseudoscalar 
resonances, X(1835) in the J/ip — >• r yrj'ir + TT~ decay process [23] and the X(pp) in the J/ip — > -fpp decay process [24]. 
It is interesting to know whether either of these has a uxf> decay mode. Based on the best solution of the PWA fit, 
new pseudoscalar states with M =1836.5MeV/c 2 , T =190.1McV/c 2 and M =1832McV/c 2 , T =76McV/c 2 are added 
in the fit, respectively. The log likelihood value improves by 2.2 and 3.5, and corresponding statistical significance is 
l.lcr and 1.6cr, respectively 

Based on the best solution, a more general test is carried out to investigate the possible contribution from additional 
resonances not listed by the PDG. Additional resonances with specified J PC and width are included (one at a time) 
in the fit, with a mass that ranges from low to high values. The scans are repeated with different widths and J PC 
values. We find that any additional state contribution has a statistical significance that is less than 5<7, and the mass, 
width and branching fraction of the A(1810) found in this way are consistent with the best solution. This method is 
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used to test whether a new resonance/state can be included in the data and no evidence for a new extra resonance is 
observed. The differences in the A(1810) parameters due to the possible presence of an additional resonance are not 
considered in the systematic error determination. 

V. SYSTEMATIC UNCERTAINTY STUDY 

For studies of the systematic uncertainties on the PWA-determined mass, width and branching fraction values for 
the X(1810), in addition to those discussed above, the effect of different background determination has also been 
studied. To estimate the systematic uncertainty associated with the background determination, the sideband regions 
(requirements III and IV) arc shifted away from the signal region by 40MeV/c 2 and 15MeV/c 2 in the ir + ir~ir° and 
K + K~ invariant masses, respectively, the side-band normalization factors are re-evaluated, and the PWA fit is redone 
using the same procedure. The differences from the best solution are taken as systematic errors. 

For the systematic errors on the branching fraction measurement, there are additional uncertainties from tracking 
efficiency, particle identification, photon detection, kinematic fit, as well as the branching fraction of the intermediate 
states and the total number of J ftp events. 



TABLE III: Summary of systematic errors 



Source 


Mass(MeV/c 2 ) 


Width(MeV/c 2 ) 


B.R.(%) 


Tracking efficiency 






8.0 


Particle Identification 






2.0 


Photon detection 






3.0 


Kinematic fit 






7.0 


Intermediate branching ratio 






1.3 


J/ip total number 






1.2 


Components in the best fit 


+2.0 
-5.1 


+12.8 
-6.7 


+19.1 
-26.1 


Resonance parameterization 


+12.2 
-1.0 


+17.3 
-33.0 


+1.6 
-40.8 


Background estimation 


+3.0 


-1.0 


-1.4 


Total 


+12.8 
-5.2 


+21.4 
-33.7 


+22.3 
-49.8 
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The systematic uncertainty associated with the tracking efficiency has been studied with J/ip — > ir + Tr~pp and 
J/rp — > KgKir, Kg — > tt + tt~ control samples [26]. The difference between data and MC is 2% per charged pion and 
kaon track. Here, 8% is taken as the systematic error for the detection efficiency of charged tracks. 

The uncertainty due to the kaon particle identification is determined from studies of a J/fi K*K control 
sample [26]. The difference in the particle identification efficiency between data and MC is 1% per kaon. Here, 2% is 
taken as systematic error for the identification of two kaons. 

The uncertainty due to photon detection efficiency is 1% per photon, which is determined from a J/ip —> control 
sample [26]. Here, 3% is taken as systematic error for the efficiency of the three photon detection. 

To estimate the uncertainty associated with kinematic fit, selected samples of ip(2S) —> ir + ir~ J/ip, J/ip — >• K + K~it° 
and ip(2S) — > it + tt~ J/ip, J/ijj — > K + K~tv°tt° events are used to study efficiency differences between data and MC. 
Compared to the final states of the studied channel, the two control samples have exactly the same charged tracks 
but one more or one less photon. The efficiency differences between data and MC are 4.2% and 7.0% for the two 
samples, respectively. Conservatively, 7.0% is taken as the systematic error associated with the kinematic fit. 

For the branching fractions of <fi — > K + K~ , ui — > 7r + 7r~7r°, and 7r° — > 77 decays, the uncertainty on these branching 
fractions listed in the PDG tables [18] arc taken as a systematic uncertainty for our measurement. The total number 
of J/ip events is (225.3 ± 2.8) x 10 6 , determined from inclusive J/ip hadron decays [9], with an uncertainty of 1.2%. 

A summary of all the uncertainties is shown in Table III. The total systematic uncertainty is obtained by summing 
up all uncertainty contributions in quadrature. The systematic uncertainties on the mass and width of the X(1810) 
arc ig 3 McV/c 2 and I34 MeV/c 2 , respectively, and the relative systematic error on the product branching fraction is 

-50 /o - 

In the best PWA fit, the threshold enhancement A(1810) is parameterized by a Breit-Wigner formula with a constant 
width. Since the enhancement structure is near the uxp threshold, other decay modes of X(1810) are expected. To 
account for this, the Flattc formula [25] is used to parameterize the structure A(1810). We assume the A(1810) mainly 
decays to ucp and K + K~ final states, and g^ and %kk are the coupling constants to the two modes, respectively. 
We test two cases, one with g^ = 1 , gxx = and the other with g^ = 0.5, gKK = 0.5, the mass and width of the 
X(1810) shift by ±19 MeV/c 2 and ±75 MeV/c 2 , respectively, and while the relative change in the product branching 
ratio is ±65.1%. These arc considered as a second systematic error due to uncertainty of the model dependence. 
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VI. SUMMARY AND DISCUSSION 



We use (225.3±2.8) x 10 6 J/tp events accumulated with the BESIII detector to study the doubly OZI suppressed 
decays of J/tp — > jujc/), uj — > ir + ir~ir°, <p — > K + K~. A strong deviation from three-body phase space for J/tp — > juj<p 
near the uxp invariant-mass threshold is observed. Assuming the enhancement is due to the influence of a resonance, the 
X(1810), a partial wave analysis with a tensor covariant amplitude determines that the spin-parity of the X(1810) 
is ++ , and the statistical significance of the A(1810) is more than 30er. The mass and width of the X(1810) 
are determined to be M = 1795 ± 7(stat)l5 3 (syst)±19(mod) MeV/c 2 and r = 95 ± 10(stat)+ 2 * (syst)±75(mod) 
MeV/c 2 and the product branching fraction is measured to be B(J/ip -> 7^(1810)) x B(X(1810) -> u<j>) = (2.00 ± 
0.08(stat)li'oo(syst)±1.30(mod))xl0" 4 , where the first error indicates the statistical error and the second is the 
systematical error. These results are consistent within errors with those from the BESII experiment [1]. 

The decay J/tp — > ^ujip is a doubly OZI suppressed process that is expected to be suppressed relative to J/tp juju 
or J/tp — > j(j)4> by at least one order of magnitude [2] . The anomalous enhancement observed at the uxp invariant-mass 
threshold and the large measured branching fractions (~l/2 of B{J/ip — > "f(f>4>) [18]) are surprising and interesting. 
The enhancement is not compatible with being due either to the X(1835) or the X(pp), due to the different mass 
and spin-parity. The interpretation of the enhancement as being due to effects of ui<p final state interactions (FSI) is 
not excluded in this analysis. Searches for this structure in different decays modes, e.g. K*K*, louj, etc., and in other 
production processes, e.g. J/tp — > <pu<p, J/tp — > Louxp etc., are essential to explore the nature of the enhancement, 
and gain more insight in the underlying dynamics. The search for other possible states decaying to uj<p would also be 
of interest. Contributions from ++ , h , 2 ++ partial waves arc found to be necessary in the PWA fit and simply 
assigned to the /o(2020), ??(2225) and /2(1950), respectively, in this analysis, since the PWA fit is not sensitive to 
those masses and widths. 
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